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Abstract 

Taking as a probe an SU(2) gauge field with Yang-Mills action in a 3 + 1 dimensional 
Lifshitz black hole background, we use the gauge/gravity correspondence to discuss fi- 
nite temperature effects in the dual theory defined on the boundary. In order to test the 
dependence of results on the anisotropic scaling exponent z we consider two analytical 
black hole solutions with z = 2 and z = 4. Apart from solving the equations of motion 
in the bulk using a numerical approach, we also apply an analytical approximation 
allowing the determination of the phase transition character, the critical exponent and 
the critical temperature behavior as a function of z. 

Models with anisotropic scaling were introduced in condensed matter physics more than 
thirty years ago in order to discuss tricritical points (see |T] and references therein). They are 
at present actively investigated in the context of gravitational theories in which space-time 
anisotropic scaling improves the short distance behavior (see |2] and references therein). A 
link between these two issues was established by Kachru, Liu and Mulligan [3] within the 
framework of the gauge/gravity correspondence by searching gravity duals of non-relativistic 
quantum field theories. Studying the equations of motion of Einstein gravity with negative 
cosmological constant coupled to p — 1 and 2-forms a solution was found in [3| with the 
metric taking the form 

ds 2 = L 2 [-r 2z dt 2 + r 2 dx 2 + ^ (1) 

where < r < oo, dx 2 = dx\ + . . . dx^, L is the radius of curvature of the geometry and 
z > 1. Metric Q is invariant under anisotropic scaling of space-time coordinates 

_ T 

t -> XH , x -> Ax , r ->• - (2) 

A 

with z playing the role of the dynamical critical exponent |2]. The coordinates' inverse 
length dimensions are: [t] = —z, [r] = +1, [x] = [y] = —1. Taking eq. ([T]) as a background 
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metric, the authors in ref. [3] extended the gauge/gravity duality to the case of models with 
anisotropic scaling and explored the boundary observables dual to free scalar fields in a 3 + 1 
dimensional bulk. 

The finite temperature extension of the gauge-gravity duality requires to consider a black- 
hole bulk metric with line element 

ds 2 = L 2 (-g z (r)r 2z dt 2 + -^—Ar 2 + r 2 {dx 2 + dy 2 )] (3) 
V 9z{r)r 2 J 

where g z vanishes at the horizon r#. Different black- hole solutions with anisotropic scaling 
are available |l]-[8] and a number of holographic studies have considered them as a back- 
ground with bulk Lagrangians including different fields: charged matter, Abelian and non 
Abelian gauge fields, and massive Proca fields [9|-[13j. 

Using the gauge/gravity correspondence we study in the present work finite temperature 
effects in the dual theory defined on the boundary We take as a probe an SU(2) gauge field 
Afj, with Yang-Mills action, this implying that the order parameter is a vector and that one 
should expect a strongly anisotropic result for conductivities (among the works cited above, 
solely ref. |13j has considered a vector order parameter). In order to test the dependence 
of results on z we shall consider two analytical 3 + 1 dimensional black hole solutions with 
different z values: the z = 2 black hole found in [8J and the one presented in |2] and |15j for 
the z = 4 case. 

The z = 2 black hole constructed in [8] arises as a solution of the equations of motion 
for a 3 + 1 dimensional gravitational theory with negative cosmological constant coupled to 
a massive vector field and a scalar field $ without kinetic term. The action reads 

S 2 = \ J d 4 x(R - 2A) - J d A x Q exp(-2$) T^ v + ^A.A^ + (exp(-2$) - 1)) (4) 

The solution of the equations of motion corresponds to a metric with line element given by 
eq. ^ with 

0»(r) = l-^f. (5) 

Starting from an action in which a Maxwell field A^ is coupled to gravity but not directly 
to the massive vector field, a charged z = 4 flat horizon black-hole solution was presented in 
refs. [9j,[15j. The action takes in this case the form 



S 4 = l -J d 4 x(R - 2A) - J d A x Qf^F^ + l -JF^ v + ^•V*") 



(6) 



with the black holes function g^ given by 



9*(r) = l-£r. (7) 



where Q an integration constant related to the Maxwell field. 
The black hole temperature associated with ^ is given by 



Tz = ~r — 4^ — (8) 
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so that for the z = 2, 4 black holes described above one has 

r 2 (l 2 
T 2 = -^, T 4 = ^- . (9) 

2vr' 8tt V 7 

Note that [T z ] = z. 

As stated above, we take as a probe an 577(2) gauge field A a ^ (a = 1,2,3) in the black 
hole background ^ with g z (r) given by ^ and ([7]). We take from here on L = 1. We start 
from the Yang-Mills action 

S = - 1 -jd i x^W\ F; v F a ^ (10) 

The field strength (a = 1, 2, 3) is defined as 

Fl v = d»Al-d v Al + e^A\Al (11) 
We have taken the gauge coupling constant equal to one. The equations of motion read 

1 d^^F a ^) - e^F^Al = 0. (12) 
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In order to solve these equations we shall consider the ansatz proposed in [16J for a relativistic 
non-Abelian gauge theory defined in an asymptotically AdS space-time 

A = <p(r)T 3 dt + w(r)r 1 dx (13) 

where r a are Pauli matrices. The gauge field inverse length dimensions are [(f)] = z and 
[a;] = 1. It will be convenient to introduce the dimensionless variable u = rn/r, so that the 
horizon is located at u — 1 and the asymptotic boundary at u — 0. In terms of this variable, 
and inserting the black hole metric ([3]), eqs.(12) reduce to 

</>" + —0' - -Y^T-M = (14) 
u rjjg{u) 

co" + ^-\d u {u l - z g{u))uj' + -^^0 2 = 0. (15) 
g{u) rfg(u) 2 

Let us discuss appropriate conditions for the gauge field components. The consistent 
conditions at the u = 1 horizon are 

0~0i(1-m) + ... (16) 

u — y 1 

U) ~ Uh + (1 — u) 2 u)i + ... (17) 

Concerning the boundary u = 0, one has 

(f) ~ lj + p\ n (u) + ... (18) 

z = 2, u — y 

CO ~ Wo 

+ ftw 2 + ... (19) 
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H + pu 2 ~ z + 



cu ~ u Q + Qu z + ... 



z > 2 , u ->■ 



(20) 
(21) 



According to the gauge/gravity correspondence p will be identified with the chemical 
potential and p with the total charge density in the dual theory defined on the boundary 
The general solution for <fi with z = 2 in the normal phase takes the form 



UJ 



Hn + p In (u) 




Using the horizon condition (f)(1) = 0, we have that 

p n = 



(22) 



(23) 



so that the chemical potential of the normal phase vanishes. In contrast, for the z 
normal phase one has 



UJ 



1 



u- 



(24) 



and hence the chemical potential of the z = 4 normal phase is non-vanishing, p n = p. 

In the z = 1 relativistic case the divergencies of the action at the boundary are eliminated 
by adding counterterms. New divergent terms arise for z > 2 but taking a fixed charge 
density p as boundary condition make these terms temperature independent jTU] . We thus 
adopt this natural choice in what follows. If, as it happens in the z = 1 case [16J, ansatz 



(13) for a z > 1 theory can be related to a an holographic p-wave superconductor, the order 



parameter should then be Q. The necessary requirement for f2 to be unsourced forces the 



choice of vanishing loq in eq. (21) or eq. (24). The divergencies of the action in the normal 
uj = phase and the superconducting w ^ one coincide leading to a finite free energy 
difference, as we shall see below. 

We shall now proceed to calculate the free energy J 7 , related to the Euclidean on-shell 
action according to 

J 7 = T z S E \ onshell (25) 



Before proceeding to the Wick rotation of the action we insert the ansatz (13) in eq. (10) 



S 



V 



2T Z 



du- 



3+z 



U 



-r 2 H Z u 2z+2 ((f)') 2 



H 



9z(u) 



(26) 



where V is the two dimensional boundary spatial volume. 

We start with the z = 2 case. Integrating by parts eq. (26) and using the equations of 
motion we get 



V 



s 



(U(j)(j)')\ u=€ - 



r 2 H g(u) 



-uj UJ\ 



u 



du 



u 



r\g{u 



J.2 2 
-0 UJ . 



(27) 
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Here e is a cut-off which will be put to zero at the end of the calculations. As discussed 
above, we choose to work in the canonical ensemble and hence we add a boundary term to 
the action [T7| 



dtd'x y /=gA fl F u "\ u= 



V 
2T 2 



u 



OJ UJ\ 



(28) 



After a Wick rotation, using eq. (25) and the boundary behavior of the gauge field the free 



energy density at fixed charge takes the form 



V 



du 



u 



r 2 H g(u) 



iV + -p 2 ln( M ) 



(29) 



The logarithmic divergent term in the r.h.s. will play no role when comparing the free 
energies of the solutions with oj ^ with that of the normal oj = case which has the same 
divergent term so that one ends with 



A J 7 T — T n 



V 



V 



-pp + 



du 



u 



r 2 H g(u 



J.2 2 

-<p OJ 



(30) 



Proceeding in the same way in the z = 4 case, we find 



V 



-p(p - p n ) + 



du 



u 



u 



J.2 2 

-<p OJ 



z = A 



(31) 



where p n is the chemical potential of the normal phase. 



Before discussing the numerical solutions of equations ( 14 )-( 15 ) we shall develop an an- 



alytic approach which allows to calculate the critical temperature and the behavior of the 
order parameter with remarkable accuracy. The method is based in a proposal presented in 
ref. [18J which consist in obtaining solutions in close form by imposing conditions of con- 
tinuity and smoothness at a point u m intermediate between the boundary (u = 0) and the 
horizon (u = 1). Originally u m was arbitrarily chosen to be 1/2 and rather good results in 
comparison with more involved numerical methods were obtained. As discussed in [19] the 
agreement stems from rather elementary considerations on perturbation of Schrodinger-like 
equations. We here extend the method in order to determine u m from a simple free energy 
argument and in this way, the method turns out to be a powerful tool to study the behavior 
of the system as a function of z. 

In practice, we shall consider expansions of the fields near u = 1 and u = and determine 
their leading orders coefficients by connecting the expansions at u — u m . We start from the 
case z = 2. For the solution near the horizon (u = 1) we have, up to second order in the 



expansions of the fields we call oj h {u) and <f> h (u), 



oj h (u) 

4> h {u) 



l 



oj^ + ul {u-l) + -0J h 2 {u-l) 2 



(32) 
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where u;f , cj)f, are constants to be determined. The superscript h indicates that the expansion 
is performed near the horizon. Now, conditions (16) at u — 1 imply that 



1 



2 + 



H 
,2 
H 



wf = , u>» 



1 <Pju) H 
16 r 4 H 



(33) 
(34) 



We now insert these relations in eq. (32) and match the expansions of oj and (f> and their 
derivatives at u — u m . From this we get 



4r| 



\A - Ur. 



U r . 



2r 



2 ^ 



1/2 



1 - Mr. 



2u m (1 "m ) 

At this point we can write r# in terms of the temperature T using eqs.([9]) 

1/2 



(35) 
(36) 



n 



i 



4ttT, 



2 - U r 



1-M m 2u m (l-w m ) 1/2 



(37) 



Determination of the point at which the order parameter Q vanishes leads to the critical 
temperature 

= _L [± !V) p , ( 38 ) 

87T M m (2 - U m ) 

One can also infer the temperature dependence of the condensate close the phase transition 

1/2 



Q = Af 2 (u m ) (4ttT 2 c 



,1/2 



T 



1 f2-U r , 



Urn \ ^™ 
2 _ Q„.4 



1/2 







Similar calculations with 2 = 4 yield, using = 8r^ 

, = Q 2 

^ (1-1^/2)1/2(1-^)1/2' " <(2-„,„) 

_ 9 i/4 ^ 4-3^ m 4p (2- Mm )V 2 \ 1/2 

r u m -i + q 2 ui(i - u m y/y 

which defines the critical temperature as 

(2- Mm ) 1 /2(l_ Mm )l/2 



(39) 
(40) 

(41) 
(42) 



T 4 C 



25/2^(4 - 3 Um ) 



-P 



(43) 
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Finally for the behaviour of the order parameter near the critical temperature we obtain 



n = M(«m) (16^ 



1/4 



M(Wm) 



1 (4 - 3u m ) 



T 

1 1 

1/2 



1/2 



(2 - Um ){i - Um y/ 2 



(44) 



(45) 



We then see that both for z = 2 and z = 4 the behavior of Q near the critical point 
reveals a typical scenario of a second order phase transition, with an ordered phase u ^ for 
T < Tl in agreement with the results in the most diverse relativistic models explored using 
the gauge/gravity duality, with critical exponents coinciding with those obtained within the 
mean field approximation, independently of the choice of u m . 

To confirm the results obtained above we have still to compare the free energy associated 
to the solution we have found with that for the disordered (normal) phase which corresponds 
to u = 0. If the difference of free energies AJ 7 , is negative below the critical temperature 
then a phase with non-vanishing order parameter will be preferred for T < T%. This fact 
will allow us to determine function of p, from minimization of AJ 7 written in terms 



of expansion (32) from the horizon to u m and of expansions (20)-(19) from the boundary to 

a t 1 r Um ii 1 r 1 a 

* -pQi-l^ + t du ^-( (j) Wy + l dw _^_(^)2 ( 46 ) 



V 



r 2 H g(u) 



r 2 H g(u) 

4. Note that we have not 
included the divergent term in ( 46 ) since we are working at fixed p and hence such term is 
u m independent. Minimization of eq. (46) gives a solution for u m which, inserted in eqs. (38) 



where and co are given by (19) for z = 2 and (21) for z 



and (43) gives the following critical temperature coefficients 



T 2 C 



0.022 p, 



Tl 



0.025 p 



(47) 



We will confirm below this scenario and compare these results with those obtained by solving 
the equations of motion numerically. Before doing this let us note that the critical tempera- 
ture obtained analytically increases when changing from the z = 2 to the z = 4 system. To 
determine whether this is a general behavior for arbitrary values of z is relevant in connection 
with the theory of Fermi liquids [TU]. To analyze this issue in more general terms one can 
take for illustrative purposes the following black hole function 



g{u;z) 



u 



(48) 



which includes the actual Lifshitz z = 2 and z = 4 black hole solutions studied previously. 
From eq. (|8|) one can write r# in terms of z and T and then, using the analytical approach 



one can confirm that, for black holes of the form (48 ), T£ is a growing function of z for z > 2 
for any choice of u m . 

We now proceed to solve the equations of motion numerically. The strategy is the fol- 
lowing: the solutions are searched as functions of the parameters Uh and <pi at the horizon 
with vanishing constant term for <fi and with general non-vanishing uq at the boundary (see 
eqs. ( 16 )-( 19 )). Then the numerical system is solved searching possible values of <fii at the 
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horizon for which Uq vanishes. In this way we have obtained a set of solutions for differ- 
ent field values at the horizon. The existence of several solutions satisfying the appropriate 
boundary conditions, each one corresponding to a different value of <pi, is a phenomenon 
already present in the relativistic case |16| . For increasing values of 0i the solution for for 
u has an increasing number of nodes n. Now, evaluation of the free energy shows that it 
increases with the number of nodes and hence we conclude that solutions with n > 1 are 
energetically disfavored so that we shall solely discuss the zero-node solution. 

Our numerical solution confirms the results found analytically: a finite temperature 
continuous symmetry breaking phase transition takes place both for z = 2 and z — 4. As 
shown in Figure 1 the system condensates at a critical temperature T c . The behavior near 
T c can be seen, by fitting the curve, to correspond to a second order transition with critical 
exponent 1/2 as advanced by the analytical result, eqs.(39)-(44). It should be stressed that 
profiles for z = 2 and z = 4 are strikingly resemblant. What distinguishes the two cases is 
the value of the critical temperatures: 

T 2 C = 0.023 p, T 4 C = 0.031 p (49) 



Comparing these values with those obtained previously using the analytic approach eqs. (47) 
we find a remarkable agreement. 




Figure 1: The numerical result for the condensate as a function of temperature for the z = 2 
(left) and z = (right) cases. The condensate goes to zero as (T - T c f/ 2 in both cases thus 



confirming the analytic results (eqs. (39), (44)) 



Note that at low temperature the condensates appear to diverge as a negative power of the 
temperature. This behavior was already encountered in the relativistic z = 1 case, both for 
s-wave |2D] and p-wave [16j holographic superconductors and can be ascribed to the relevance 
of back-reaction when the condensate becomes too large so that the probe approximation is 



no more valid. Using again eq.(48) as an illustration, our analytical approach shows that the 



behavior of the condensate for T small -in the range of validity of the probe approximation- 
is f2 oc T~( z ~ 2 ^ 2z for z > 2 independently of the choice of the matching point u m . 



Using formulae (30)- (31 ) we have computed numerically the free energy difference between 
the ordered and disordered phases (see Figure 2) confirming that, both for z = 2 and z — 4, 
the ordered phase is preferred below the critical temperature whose values coincide with 



those given by (49). 
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Figure 2: The free energy difference between the condensed and the uncondensed phase as 
a function of temperature for the z = 2 (left) and the z = 4 (right) models 



Finally, we shall compute the electromagnetic response to small time dependent pertur- 
bations of the electromagnetic fields in the ordered phase. To do this, we start from the gauge 
field ansatz (13) (that we shall denote A°J d (u) for clarity) and following [16j we consider the 
perturbation 



A, 



A° rd (u) + a,(u,t) 



t d x v = e -*»f* [(a 1 ^ 1 + a 2 t T 2 )dt + a 3 x T 3 dx + a 3 y T 3 dy] 



(50) 
(51) 



with wj the frequency associated to the perturbation. The linearized Yang-Mills equations 
read 

ia^P") - e abc F^A^ - t abc F^a, c = 



-g 



where 



^.abc 



Using eqs. (13), (51) one finds four second order equations 
1 — z g'{u 



3" , / tz^. , \ 3 ' , w ) u2z 2 a s _ 



u 



9W 
z , g'(u) 



co- 



ll 



9W 



u 2z-2 



a a = o, 

r 2 H g(u) y 



a",. ( — | '/',' + Q , n/ N (-iwfua 2 + wja 3 x - uupa]) = 0, 



a' + 



Z ~K i 



2" , 

a; + 



u 

z-1 



at + 



U(j) 



r 2 H g(u) 



to 



u 



9 / \ t 

rjjg[u) 



r 2 £g 2 {u) 



ZWfU} -a 3 = 0. 



r 2 H g(u) x 



and two first order equations 



iwfd] + (f)a 2 — 4>a 2 = 
iwfaf — (pa]' + (f)'a\ — g(u)u 2 ~ 2z (ujd u — u)')a 3 = 0. 



(52) 
(53) 

(54) 
(55) 
(56) 
(57) 

(58) 
(59) 
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Let us concentrate on the case z — 2. The choice of the electromagnetic perturbation should 
correspond to a wave traveling away from the conformal boundary at u = (an "in-going" 
wave). In the present case one has, near the horizon 



iw j? ■ w f 

^ — ^Yl *,2\ 2r^ (a | \ 3 _ D(A n .\ Sri 



a 3 = a{\ -v?) 2r H {! + ...) , < = /3(l-u) 2r « (1 + a x {l - u) + . . .) (60) 

a\ = 7 (1 - (a 2 (l - w) 2 + . . .) , a 2 = 5(1 - u)~**t (a 4 (l -«) + ...) (61) 

with a, . . . , 7 dimensionful constants. At the boundary we have instead 

a 3 y = a 3 m + u 2 al {2) + ..., a 3 x = a 3 x{0) + u 2 a 3 x{1) + ... (62) 

4 = a l(o) + «f(i) ln («) + ••■> a 2 = a 2 (0) + a 2 (1) ln(w) + ... (63) 

where all coefficients <2j can be determined as functions of a; and at the horizon once Wf is 
specified. 

The conductivity can then be obtained using Ohm's law. Following [16J for the case of 
non-Abelian gauge fields, the conductivity components are 

a „ = -i*%L, a _ = -J%-U )+Q *»*f> + '? '<»> ). (64) 

We show the numerical solution for the real and imaginary parts of a xx and o yy for the 
z = 2 system in figures 3. As in the relativistic case the conductivity components approach 
1 at large Wf. We observe the formation of a gap in the real part of a yy as it happens 
in the case of a Maxwell field coupled to a scalar [20] and in the purely Yang-Mills |16| 
bulk Lagrangians cases. There is a pole in the imaginary parts of a xx and a yy at Wf = 
characteristic of superconducting behavior. There is a second pole in the imaginary a xx at 
Wf = = 0.199p at T/p = 0.022 accompanied by the corresponding delta function in its 
real part, in agreement with Kramers-Kronig relations (this delta function is not represented 
in figure 3 left since the numerical procedure can only render continuous functions). The 
value obtained numerically satisfies = fi as expected from eq.(64). This pole is absent in 



the analysis of |13j for a bulk Yang-Mills Lagrangian in the background of a different z = 2 
Lifshitz black hole (the one presented in [3] with Q2{u) = (1 — u 4 ) arising in the case in which 
the dilaton field is dynamical, instead of the one we have used, eq.([5])). In [13] such absence 
was attributed to the logarithmic behavior of Aq resulting from the z = 2 scaling. Our result 
shows that for the z = 2 black hole background that we used such logarithmic behavior does 
not prevent the existence of this pole. 

The analysis of the z = 4 theory follows similarly and the behavior of conductivity 
components is qualitatively the same. We also find in this z = 4 case, with g^{u) = 1 — w 4 , 
a second pole located at = 20. 5p for T/p = 0.022. 

We shall end this work with a brief summary and a discussion of our results. We have 
studied finite temperature effects in two models with different dynamical critical exponent 
using the gauge/gravity correspondence. Looking for a vector order parameter and inspired 
by Gubser and Pufu's work on z = 1 p-wave holographic superconductors p.6j, we have 
chosen as gravity dual a Yang-Mills theory in the gravitational background of Lifshitz black- 
holes with z = 2 and z = 4. Apart from solving the equations of motion in the bulk 
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using a numerical approach, we have also extended the analytical approximation developed 
in [18j-[19j which allows to reproduce the numerical results with remarkable simplicity and 
precision. 




Figure 3: Real an imaginary parts of conductivity as a function of the frequency for T/p = 
0.022 for the z = 2 system. The solid line corresponds to a xx and the dashed one to a yy . 
The insert figure on the right displays a detail of the imaginary part of a yy rendering visible 
the pole at cjf = 0. 

Although one could presume that the anisotropic scaling of the background metric would 
lead to a critical behavior differing from the one found in [16] for z = 1, our results show 
instead a remarkable resemblance with the relativistic case. In particular, the condensate 
has the typical (T£ — T) 1 / 2 mean field behavior for T close to the critical temperature T£ 
both for z = 2 and z — 4. The dependence on z only affects the coefficient in the critical 
temperature which grows with z, a behavior that could be argued to be valid for arbitrary z, 
as we illustrated applying our analytic approach to an heuristic black hole function g(u; z) 



defined in eq.(48). Using the same approach we were able to extract the condensate behavior 
in the range of small temperatures where the probe approximation is valid, finding that 
Vt oc T _(z ~ 2 )/ 2z , in total agreement with the numerical calculations. All these results confirm 
that the analytic approximation developed in [18j-[19j and refined here has proved to be 
sufficiently accurate as to avoid the necessity to resort to numerical methods. 

Acknowledgments : This work was supported by CONICET , ANPCYT , CIC and UNLP, 
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